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Abstract 

We identify a class of perturbatively computable measures of interjet energy flow, which can be 
\ associated with well-defined color flow at short distances. As an illustration, we calculate correlations 

between event shapes and the flow of energy, Qq, into an interjet angular region, fl, in high-energy 
two-jet e"'"e~-annihilation events. Laplace transforms with respect to the event shapes suppress states 
' with radiation at intermediate energy scales, so that we may compute systematically logarithms of 

CNJ I interjet energy flow. This method provides a set of predictions on energy radiated between jets, as a 

function of event shape and of the choice of the region SI in which the energy is measured. Non-global 
logarithms appear as corrections. We apply our method to a continuous class of event shapes. 
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P-i. 1 Introduction 

O , Inclusive jet cross sections have a long history, and the agreement of predictions based on perturbative 
QCD with experiment is often impressive. The study of energy flow into the regions between energetic 
jets [1] gives important, and in some sense complementary information about the formation of final 
^ . states. Interjet energy flow is expected to reflect the flow of color at short distances [1, 2, 3] and 
^ I may give insight into the hadronization process. In addition, knowledge of the relation between energy 
and color flows is important for investigating hard collisions [4], to distinguish QCD bremsstrahlung 
originating in the underlying hard scattering or decay from radiation induced by multiple scatterings. 
The study of observables in jet events that test the flow of energy is thus of considerable interest. The 
computation of such observables, however, has turned out to be subtle [5], for reasons that we will 
briefly review below. For deflniteness, we will discuss these issues in dijet events at large center-of-mass 
energy Q in e"^e~ -annihilation [6]. In this case the underlying color flow at short distances is unique, 
corresponding to the creation of a quark-antiquark pair. This enables us to illustrate our method, 
while avoiding the complication of multiple color exchanges, for example, singlet and octet exchange in 
quark-antiquark scattering [2, 3]. 

We study the flow of energy or transverse energy, Qq = eQ, with e <C 1, into an interjet angular 
region, $7. The process 

e+ + e~^ 2 Jets +X^ + Rn{Qn), (1) 

'^Based on talk given by C. F. Berger at TH-2002, International Conference on Theoretical Physics, Theme 2: "QCD, 
Hadron dynamics, etc.", Paris, France, 2002. 
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is illustrated in Fig. 1. represents soft radiation into the region between and the jet axes, denoted 
by SI, and denotes soft radiation into S7. Our goal is to relate soft radiation into region Vt to the 
hard scattering that produces the jets. 




Figure 1: Sources of global and non-global logarithms in dijet events. Configuration 1, a primary 
emission, is the source of global logarithms, configuration 2 results in non-global logarithms. 

Following [5], we refer to observables as "non-global" when they are defined by radiation into a spe- 
cific portion of phase space, while remaining more-or-less inclusive elsewhere. In non-global cross sections 
there are two sources of logarithmic corrections, as indicated in Fig. 1: "primary" emissions, such as 
gluon 1 in Fig. 1, are emitted directly from the hard partons into Vl. Phase space integrals for these 
emissions contribute single logarithms per loop: a" ln"((5/(5n) = a"ln"(l/e). These logarithms expo- 
nentiate, and may be resummed in a straightforward fashion [3]. There are also "secondary" emissions 
originating from the complementary region 0, illustrated by configuration 2 in Fig. 1. As emphasized 
by Dasgupta and Salam [5], emissions into from such secondary partons can give logarithms of the 
form ln"'(QQ/QQ), where is the maximum energy of radiation in Ct. These have become known 
as non-global logarithms. If no restriction is placed on the radiation into Cl, then can approach Q, 
and the non-global, secondary logarithms can become as important as the primary logarithms. The 
non-global logarithms arise because real and virtual enhancements associated with secondary emissions 
do not cancel each other fully at fixed Qq,. In general, such non-global logarithms, as studied in [5] and 
[7], are sensitive to color flow at all scales and in all directions, and they do not exponentiate in the 
same simple manner as do logarithmic enhancements from primary emission. In a sense, they mask the 
underlying color exchange at short distances. 

In the studies of non-global observables in Refs. [5, 7] , numerical estimates of both types of logarithms 
were carried out in various situations. Here we will adopt a somewhat different strategy, and introduce 
correlations between energy flow and event shapes ("flow/shape correlations"). Our aim is to select final 
states whose formation is sensitive primarily to radiation from the highest-energy jets. By restricting 
the range of event shapes, we can systematically limit radiation in region while retaining the chosen 
jet structure. This allows us to treat all leading effects analytically. In essence, the color and energy flow 
originating from the primary radiation is less obstructed by secondary radiation. Our approach does 
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not resum non-global logarithms, but allows us to address the question of interjet energy flow within a 
more global context, while retaining the original motivation and physical picture described above. 



2 Energy Flow/Event Shape Correlations 

The following observation allows us to avoid leading contributions from non-global conflgurations such 
as gluon 2 in Fig. 1. Large non-global logarithms originate from emissions into Q at relatively large 
angles from the jets, because radiation from partons close to the jet directions is emitted coherently. 
By fixing the value of an event shape near the two-jet limit, we avoid final states with large energies 
in Q at fixed angles from the jet axes. To impose this two-jet condition at fixed energy flow, we will 
introduce shape functions fc, c = 1, 2, one for each jet, for radiation into D,. The relevant cross section 
for the process (1), at fixed direction for jet 1, hi, is then defined as 

1^ = ^LlM(iV)P(2.)nV-») 

X S{e - f(N))S{i - /i(JV) - /2(JV)) i(r>i - n(JV)) , (2) 



with the energy flow fixed by 



We sum over all final states A'' that contribute to the weighted event, with M{N) the amplitude for 
e+e~ N. The total momentum is pi, with pj = Q^. To define the total event shape e in Eq. 
(2) wc divide Q into two hemispheres, each surrounding one of the jets. Each hemisphere contributes 
separately, through fi and /2. An example of a suitable shape function is the thrust-like weight 

MN) = ^ E ^ii^-cosOi), c=l,2, (4) 

where 9i are the angles with respect to the jet directions. The cross section Eq. (2) measures the 
correlation oi e = fi + f2 with the energy flow into fi. Since we are interested in two-jet cross sections, 
we flx the constants e and e to be much less than unity: 

<£,£<!. (5) 

In this limit, we can neglect recoil due to the soft radiation. Generalizations to more than two jets in 
the final state, and to correlations with transverse energy flow into 17, are clearly possible. 

In the limit (5) the cross section (2) has large logarithmic enhancements in Inl/e and Inl/e, which 
we will resum. However, we will not resum logarithms like ln(e/£). Instead, the following generalization 
of Eq. (4) to the weights fc{N, a) allows us to study correlations of jet structure with energy flow into 

n-. 

fciN, a) = ^ J2 sin" Oi (1 - cos , c = 1, 2, a < 2. (6) 

As a — >■ 2 the weight vanishes only very slowly for 0i 0, and at fixed fc, the jet becomes very 
narrow. On the other hand, as a ^ —oo, the event shape vanishes in nearly every direction, and the 
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cross section at fixed fc becomes more and more inclusive in radiation into Cl. In this limit, non-global 
logarithms reemerge as leading effects. For a = we obtain the thrustlike weight (4) discussed above; 
the case a = 1 corresponds to the jet broadening. The effect of a on the shape of the radiation into i). 
is illustrated in Fig. 2. In Fig. 2 we compare the phase space available to a particle at fixed e for three 
different values of o. The radial magnitude of each plot is the maximum energy uj found from Eq. (6): 
r = eQ sm~"' 0(1 — cos0)"~^. For a = 1, as shown in Fig. 2 a), the particle is restricted to be close to 
the jet axes, while for a = 0, and a = —1, depicted in Figs. 2 b) and c), respectively, the particle is 
allowed to be farther away from the axes. For a > 1, recoil of the jets against soft radiation cannot be 
neglected [8], so that our analysis applies in the range a < 1. With the above choices of weight functions 
the cross section (2) is infrared safe, and logarithms of e and e may be resummed. 




Figure 2: Illustration of the effect of the parameter a in the weight (6) on the shape of the event: a) 
shape for a = 1, b) shape for a = 0, c) shape for a = —1. The jet axes are in the vertical direction. The 
radial normalization {eQ) is arbitrary, but the same for all three plots. 



3 Factorization and Resummation 
3.1 Factorization 

Following the procedure developed and extensively discussed in Refs. [9, 10, 11, 12], we first identify the 
"leading regions", which can give rise to logarithmic enhancements associated with lines approaching 
the mass shell. When we sum over all diagrams that have a given final state, the contributions from 
these leading regions may be factorized into a set of functions, each of which corresponds to one of these 
generic subdiagrams. The cross section, Eq. (2), then becomes a convolution in e, between functions 
Jc that represent the evolution of the primary partons, and a function S that represents coherent soft 
radiation: 

dededfii ^ ^ ?i, 6, "-i, /«) J dss S{e,es,a, fi) 

xJJ jdej,Jc{ej,,pj,,ic,a,n)5{e-ej^-ej2-£s), (7) 

c=l 
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with pj^ the jet momenta. jJL is the factorization scale, which for simphcity we set equal to the renormal- 
ization scale, dao / dfi\ is the Born amplitude for the production of a single particle (quark or antiquark) 
in direction fii. The 'short-distance' function H(Q, S,i,C2,fi'i, fJ') = 1 + C'(as(^^)), which describes correc- 
tions to the hard scattering, is an expansion in ag with finite coefficients. In Eq. (7) we have suppressed 
some arguments of the soft function that are not relevant for the discussion below. The factorized form 
(7) is illustrated in Fig. 3. 

In (7), the "jet functions" Jc model radiation collinear to the primary partons, with directions +ni 
for jet 1, and —hi for jet 2, respectively, and are constructed to be independent of e, which measures the 
radiation at wide angles into il. The vectors are introduced in the factorization of the jet functions 
from the hard scattering with the help of Ward identities. For further details we refer to [12] and to [6], 
where the jet functions arc related to QCD matrix elements, following [11]. 

The soft function S describes soft radiation at wide angles from the jets, into 0, and 0, alike, and 
thus depends on both e and e. From the arguments given in [3, 12], it follows that radiation at wide 
angles from the primary hard partons decouples from the jets, and can be approximated by an eikonal 
cross section a^'''^^ built out of path-ordered exponentials 

$«(oo,0;x) ^ p^-^9jrd^f^-^''H>^P+^)^ (8) 

where the /3's are light-like velocities in the directions of the jets. The superscript (f) indicates that 
the vector potential takes values in representation f, in our case the representation of a quark or an 
antiquark. An eikonal cross section also contains enhancements for configurations collinear to the jets, 
which in the factorization (7) are already taken into account in the partonic jet functions. Therefore, the 
eikonal cross section itself is not a suitable soft function. Rather, to avoid overcounting and to include 
only soft, but not collinear enhancements, the soft function is defined through the refactorization [2] 

a^"'^^ (e, £eik, a,fi) = j dig S (e, ig, a, ji) 

X n / ^'""^ (£c, a, n) 5 (£eik -Ss-ei- £2) ■ (9) 

c=l 

In Eq. (9) we have refactorized the eikonal cross section in the same manner as the partonic cross 
section (7), into eikonal analogs of jet functions, and the same soft function [6]. We have omitted the 
dependence of the functions in (9) on the eikonal lines f5c and for better readability. 



3.2 Resummation 

To disentangle the convolution in (7), we take Laplace moments with respect to e. 



de dhi Jo de de dhi 

= ^ i/(6,&,ni,g) n -^^(^^'^'"'"^^^O S{e,u,a,as{Q)) , (10) 

where we have exhibited the relevant ^c-dependence of Jc [6]. Here and below unbarred quantities are the 
transforms in i, and barred quantities are untransformed. In Eq. (10) we have set the factorization scale 
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Figure 3: Factorized cross section (7). The vertical line denotes the final state, separating the amplitude 
(to the left) and the complex conjugate amplitude (to the right). 



to Q, which avoids large logarithms in the hard function. Large logarithms of e and or e then occur in 
the soft and jet functions. In addition, the jets contain potentially large logarithms of pj^ • £,c/ (/^ \ic\)- 

The resummation of these logarithms follows from the independence of the physical correlations 
from the factorization scale and the choice of eikonal factors 



d 



d da (e, v, a) 
dn de dhi 
da (e, u, a) 



0, 
0. 



(11) 
(12) 



Following the methods described in Refs. [11, 12, 13], we can use (11) and (12) to derive evolution 
equations for the jet and soft functions. From this procedure, we derive a resummed expression for the 
Laplace-transformed correlations [6], which may be written as 



da (e, u, a) 



dsdh 



^ i7(^i,6,ni,Q)5 (ez^, a,0,as(eQ)) exp 
ani 



Q/2 

/ y7.(«.(A)) 



Q/2 



yTJc MX)) 



Q/2 



X exp 



d\ 
T 



2i.l/(2-a) 



2,.l/{2-a) 
C2A 



5c (a, a, (A)) + j 



dX' 

—Ac (a,a,(A')) 



Ci 



(13) 



where we have exhibited the relevant arguments of the soft function. Convenient choices for the scales 
in the A' integral are C2 = 1 and Ci = exp[— 7^; + (1 — a)/2]. The functions Ac and Be are associated 
with the jet functions [11, 13], and the anomalous dimensions 7^ and 7j^ belong to the soft and jet 
functions, respectively [6, 11]. At lowest order, A^ = 2Cc{as/iT), with Cc = Cp for quarks and Ca for 
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gluons. Other explict forms are given in [6]. Only the soft function retains dependence on the geometry 
of the region Q, as exphcitly indicated in (13). Corrections in Eq. (13) of the form as(Q) In 1/(61^) 
occur entirely within the soft function S. For large enough and/or [a|, at fixed e, large non-global 
logarithms reemerge in that function, which may be computed in the eikonal approximation. 

So long as we are able to compute corrections in as{Q) In l/(ez/) perturbatively, we are free to choose 
between eQ and Q/v in any term in (13) that generates single logarithms. Of particular relevance is 
the 7s integral in the first exponential of (13), which reflects the choice of renormalization scale in the 
soft function S. At higher orders in S, logarithms will arise from gluon emission into angular regions for 
which the upper limit on energy flow does not match the renormalization scale. Within Q the upper limit 
is eQ, and outside, in Q, it is Q/i'. (Recall, there are no collinear enhancements in S.) These logarithms 
appear multiplied by the size of the angular region from which they arise. When we measure energy flow 
into a region f2 of small angular extent, it is more convenient to choose a renormalization scale Q/i^, to 
avoid generating logarithms of eiy multiplied by the size of the larger region 0. Correspondingly, when 
subtends the bulk of the unit sphere, it is convenient to make the choice of eQ, as in Eq. (13). 



4 Results 

Eq. (13) resums single logarithms of e and single and double logarithms of v. The latter stem from the 
double integral in the last exponent on the right-hand side. In Figs. 4 and 5 we show some examples of 
numerical results from (13). 




Figure 4: Differential cross section ^ dao/dnT ' normalized by the lowest order cross section, at Q = 100 
GeV, as a function of e and a at fixed u: a) u = 10, b) = 50. is a slice (ring) centered around the 
jets, with a width of Ary = 2. 



Fig. 4 shows the dependence of the differential cross section (13), multiplied by e and normalized 
by the Born cross section, ^^^"j^^^i on the measured energy e and on the parameter a, at fixed v. 
We choose the region to be a ring around the jets, centered in their center-of-mass, with a width of 
Ar/ = 2. The center-of-mass energy Q is chosen to be 100 GeV. In Fig. 4 a), we plot ^'^^l^fdhi ^ ^^'^ 
v = 10, in Fig. 4 b) for v = 50. The cross section falls as a increases, where the jets are restricted to be 
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very narrow. Similarly, as v increases, the radiation into the complementary region is more restricted, 
as illustrated by the comparison of Figs. 4 a) and b). On the other hand, as |a| increases at fixed e, the 
correlations (f3) approach a constant value. For a — > — oo, however, non-global dependence on e and 
\a\ will emerge from higher order corrections in the soft function. 

To illustrate the sensitivity of these results to the flavor of the primary partons, we study the 
corresponding ratio of the correlation to the cross section for gluon jets produced by a hypothetical color 
singlet source. Fig. 5 displays the ratio of the differential cross section da'^{£, a) / {dedfii) , normalized by 
the lowest order cross section, to the analogous quantity with gluons as primary partons in the outgoing 
jets, again at Q = 100 GeV. We multiply the ratio by Ca/Cf to compensate for the difference in the 
normalization of the lowest order soft functions. Gluon jets are wider, and hence are suppressed relative 
to quark jets as a or increases. Figs. 5 a) and b) exhibit this behavior, where we compute the cross 
section for the same ring around the jets, centered in their center-of-mass, with a width of Ar/ = 2. 
Again, Fig. 5 a) is at = 10 and b) is at = 50. These results suggest sensitivity to the more complex 
color and flavor flow characteristic of hadronic scattering. 




Figure 5: Ratios of differential cross sections for quark to gluon jets ^ ( ^^'^1'^/!^"^^ ) ( ^'^'dJ^^fdih^ ^ ) 
at Q = 100 GeV as a function of e and a at fixed ly: a) u = 10, h) v = 50. fi, as in Fig. 4, is a slice 
(ring) centered around the jets, with a width of Ary = 2. 



5 Conclusions 

We have described a study of interjet energy flow/event shape correlations in e^e'^ dijet events, where 
all leading effects are treated analytically. By weighting the final state radiation appropriately, our for- 
malism is sensitive mainly to radiation stemming directly from the primary hard scattering. Transforms 
in a class of weight functions enable us to control the influence of secondary, or non-global, radiation on 
the energy flow between the jets, and in principle to study correlations between jet structure and energy 
flow. The study of next-to-leading order and higher corrections in which these effects appear will also 
be of interest. The application of our formalism to multijet events and to scattering with initial state 
hadrons is certainly possible, and should shed light on the relationship between color and energy flow 
in hard scattering processes. 
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